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ABSTRACT 

With the steady increase of the sample size of observed microlenses towards the central 
regions of the Galaxy, the main source of the uncertainty in the lens mass will shift from 
the simple Poisson noise to the intrinsic non-uniqueness of our dynamical models of 
the inner Galaxy, particularly, the Galactic bar. We use a set of simple self-consistent 
bar models to investigate how the microlensing event rate varies as a function of 
axis ratio, bar angle and velocity distribution. The non-uniqueness of the velocity 
distribution of the bar model adds a significant uncertainty (by about a factor of 
1.5) to any prediction of the lens mass. Kinematic data and self-consistent models 
are critical to lift the non-uniqueness. We discuss the implications of these results for 
the interpretation of microlensing observations of the Galactic bulge. In particular we 
show that Freeman bar models scaled to the mass of the Galactic bulge/bar imply a 
typical lens mass of around 0.8Mq, a factor of 3-5 times larger than the value from 
other models. 
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1 INTRODUCTION 

Five years of searches for gravitational microlenses, which 
significantly increase the brightness of a background source 
star when they come to the same line of sight by chance, have 
produced roughly 200 possible microlensing events towards 
the Galactic bulge/bar from several surveys (e.g., Udalski 
et al. 1994; Alard et al. 1995; Alcock et al. 1997; The ex- 
pectation is that with more and more events coming in at 
a steady rate, we will obtain a better understanding of the 
structure of the Galaxy, and of the mass spectrum of lenses 
(see Pacyznski 1996 for a review). Microlensing is perhaps 
the only direct method to probe the mass function in the 
Galactic bulge/bar, in particular the fraction of low lumi- 
nosity objects (brown dwarfs/M-dwarfs) just below or just 
above the hydrogen burning limit. These objects are so faint 
intrinsically {My > lOmag) and so far away (distance mod- 
ulus of about 14.5 mag) that they are missing in even the 
deepest star count studies of the bulge/bar with the Hubble 
Space Telescope (Gould, Bahcall & Flynn 1996). 

Several attempts have been made to estimate the typ- 
ical mass of the observed microlenses (Kiraga & Paczyriski 
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1994; Zhao, Spergel & Rich 1994; Han & Gould 1996; Mao 
& Paczynski 1996) on the basis that events last longer for 
more massive lenses, with m cc t^, where m is the mass of a 
single point lens and t is the time for a source on the Ein- 
stein ring (where the amplification factor is 1.34) to move to 
exactly behind the lens. The time scale t can be derived by 
fitting the light curve of the amplified source, but to convert 
it to the lens mass m requires knowing other system pa- 
rameters, including the distances and transverse velocities 
of the lens and the source. These, unfortunately, are often 
poorly known from observations. Nevertheless the problem 
can be partially circumvented if one has a well-determined 
dynamical model of the inner Galaxy. The missing informa- 
tion can then be simulated by drawing random samples of 
the lens and the source in a Monte-Carlo fashion from the 
model phase-space distribution of the inner Galaxy. 

With the increase of the event sample size, it becomes 
meaningful to ask the question whether the whole mass spec- 
trum of the lenses can be determined in the above way (Han 
& Lee 1997). This requires an understanding of the uncer- 
tainty of the underlying dynamical model. Models of the 
inner Galaxy are, unfortunately, still far from being well- 
determined and they are subject to constant modifications 
and improvements, driven by new observations (see e.g., re- 
views by de Zeeuw 1993; Gerhard 1996). How much the 
Galactic bar differs from a simple oblate rotator is an un- 
settled issue. Generally speaking the density distribution of 
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Figure 1. Two simple discs. Models to the right (R-type model) 
are rotating models where all stars (lenses and sources) move in 
the clockwise direction, and models to the left are non-rotating 
with half of the stars (lenses and sources) orbiting in the clock- 
wise sense, half anti-clockwise. The transverse velocity of the lens 
and the source at the line of sight to the center (dotted lines) is 
iiVc- The lower panel draws schematic distributions of the rela- 
tive source-lens speed V = \vs — vi\ in the simplest case where 
the lens and the source are at the front or on the back side of an 
annulus. 

the Galactic bar is constrained only up to a one-parameter 
sequence by the integrated COBE map, with the bar angle 
being a free parameter (Binney, Gerhard & Spergel 1997; 
Zhao 1997). Various velocity structures of the bar could also 
be consistent with the same bar potential (Pfenniger 1984). 
Schwarzschild-type models for the inner Galaxy (Zhao 1996; 
Binney et al. 1997) which attempt to fit simultaneously the 
photometric and kinematic data measurements are still rare. 

The observed velocity field of a bar is the result of its 
tumbling motion, the internal streaming motion on top of 
the pattern rotation, and the velocity anisotropy. Bars of 
different pattern speed and/or orbital compositions are in- 
distinguishable in their projected density and optical depth 
maps, but they differ at the level of event rate and event du- 
ration distribution. This is easily seen from the basic equa- 
tion of microlensing, 

l-l-H'^r"' a. 

where G is the gravitational constant, c is the velocity of 
light. Re is the Einstein ring radius, t is the event duration 
(one Einstein radius crossing time), Di and Ds are the dis- 
tances to the lens and the source, V is the relative transverse 
speed between the two, and m is the mass of the lens. For 
most microlensing events one observes only their duration. 



t, from which one can obtain a lensing probability (optical 
depth) T and an event rate F observationally. Other sys- 
tem parameters, such as the speed V and distances Ds and 
Di, can at best be inferred statistically from a phase-space 
density model of the Galaxy. As a result the rate is related 
to the dynamical model as follows (cf. eq. [7] of Paczynski 
1986) 

where the brackets indicate averages. In order to emphasize 
the effect of the velocity distribution we have fixed the mass 
of the lenses m and the positions of the lens and source. 

Eq. (^) can also be interpreted as a constraint on lens 
mass. If the event rate is fixed at the observed value and the 
lens volume density model is fixed to fit the observed optical 
depth, then eq. (El) implies that 

moc{Vf, (3) 

so that the inferred lens mass is sensitive to the velocity 
distribution of the dynamical model. An uncertainty of 1.2 — 
2 in velocity translates to a factor of 1.4 — 4 in the lens mass, 
in which case one could mistake, e.g., an M-dwarf dominated 
model for a brown-dwarf dominated model. 

To illustrate this point, it is best to take one step back- 
wards to consider the simplest axisymmetrical models. If one 
observes the two discs shown in Figure |l] edge-on, then the 
distribution of transverse relative source-lens speed peaks at 
V = 2Vc for the maximum rotating disc (R-type; R stands 
for rotator), but peaks at both V — 2Vc and V = Q for 
the two-stream counter rotating disc (NR-type, standing for 
non-rotator). Since the median speed (y) is twice as large in 
the R-type model as in the NR-type model, twice as many 
events with duration twice as short are expected in the R- 
type model as in the NR-type model. In a more realistic 
situation where lenses and sources are distributed across a 
uniform disc and they orbit the center at constant angular 
speed, the R-type model would predict (after a simple cal- 
culation) an average speed (y) and an event rate F a factor 
of 1.2 times longer than in NR-type models. 

Past attempts to statistically determine the mass func- 
tion of microlenses all suffer from uncertainties in the dy- 
namical models of the inner Galaxy. Binney, Gerhard & 
Spergel (1997) and Zhao (1997) studied volume density 
models of the Galactic bar consistent with the projected 
light distribution from the COBE/DIRBE maps, Zhao & 
Mao (1996, hereafter ZM) addressed the effects of the un- 
certain volume density model on the microlensing probabil- 
ity (the optical depth). This paper extends such studies to 
the microlensing event rate distribution. We concentrate on 
breaking the degeneracy of dynamical models with the same 
three-dimensional density distribution but different velocity 
fields, which are indistinguishable in studies based on the 
projected light and the optical depth. 

The main goal of this paper is to study the trend of the 
microlensing event rate as a function of the key parameters 
of the bar (the pattern speed, the axis ratio and the orien- 
tation angle) assuming a fixed lens mass m. This helps us 
to gauge whether the lens mass should increase or decrease 
in case a modification of the bar parameters is required to 
be consistent with some future observations or dynamical 
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models. A side result is to give an error bar for predictions 
of the lens mass due to the non-uniqueness of the bar. 

The rigorous way to survey microlensing properties of 
Galactic bar models is to build many sequences of three- 
dimensional bar models such as in Zhao (1996) which cover 
the multi-parameter space. But such heavy numerical mod- 
elling would be very inefhcient. For a first study it is more 
interesting to gauge the underlying physical effects with sim- 
ple theoretical models. In this paper we show which insights 
can be gained from studying the microlensing properties of 
the analytical two-dimensional Freeman (1966) bars, which 
surprisingly capture most of the microlensing effects of a 
realistic bar fairly well. 

The outline of the paper is as follows. In §2, we sum- 
marize the properties of the Freeman bars, and we derive 
the optical depth for self-lensing. In §3, we derive basic scal- 
ing relations for the event rate distribution. In §4 we show 
the run of event rate as functions of bar parameters. In §5 
we show how to distinguish Freeman bars with the same 
projected density and optical depth map, and in §6 we de- 
rive the rate for extremely short events. Finally, we discuss 
implications for the Galactic bar in §7. Some mathematical 
details are given in the Appendix. 



2 SELF-LENSING OF THE FREEMAN BARS 

Freeman (1966, hereafter F66) discovered self-consistent 
tumbling bar models with a known analytical distribution 
function. These two-dimensional bars have been widely used 
to gain insight into the structure of general self-consistent 
bars (Hunter 1974; Tremaine 1976; Weinberg & Tremaine 
1983). Despite the two-dimensional nature and the special 
distribution functions, this class of easy-to-build bars con- 
tains the main factors by which bars can affect the mi- 
crolensing rate. These include its elongated density contours 
and similarly elongated velocity ellipsoid, plus the pattern 
rotation and streaming motion. Our aim is to vary these fac- 
tors and "observe" the bar both at a variable line-of-sight 
angle from its long axis and at a variable projected radius. 



2.1 Freeman bars 

The Freeman bar is a two-dimensional inhomogeneous ellip- 
tical bar of finite extent, with surface density E given by 



and a potential 
GM 



a2 &2 ' 



2{ahy- 



■[A2{q)X' + B2{q)Y^ 



1= 



(4) 



(5) 



where M is the total mass, a and 6 are the semi-major (X) 
and semi-minor (Y) axes of the bar, and A2{q), B2{q) are 
dimensionless function of the axis ratio q given in equation 
(2) of F66. The symmetry axes {X and Y) rotate with an 
angular speed Q with respect to the rest frame. For a given 
density, there is a sequence of self-consistent models with dif- 
ferent amounts of pattern rotation. The sequence is a func- 
tion of fl/Qj, where 

''A2- 
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Figure 2. Orbits in Freeman bars with GM = 1 and semi-axes 
a = 1.420 and b = 1.084 and a boundary as indicated with the 
dotted ellipses in each panel. The orbits are drawn in the X — Y 
frame corotating with the tumbling bar with a pattern rotation 
speed Q = 0.07 (left panels) and 0.67 (right panels). All orbits are 
assigned the same Jacobi integral Ej = 0.5, which confines the 
orbits within the zero velocity curves shown as the dashed eUipses 
in each panel. The heavy solid lines indicate the boundaries of 
each orbit, as given in Appendix A. 



is so defined that when Q = flj the bar is of Jacobi type, 
that is, with neither internal streaming motion nor velocity 
anisotropy. 

The potential (P) is that of a two-dimensional 
anisotropic harmonic oscillator with generally incommensu- 
rable frequencies in the X and Y directions. The motion is 
separable in canonical coordinates {Pi, P2,Qi,Q2) defined 
by F66. Accordingly, all orbits have two independent iso- 
lating integrals of motion, and they are in fact rectangular 
Lissajous figures in the (Qi, (32)-plane. In the corotating 
{X, Y')-frame, the orbits show a more interesting variety of 
shapes, similar to what is seen in more general bar poten- 
tials. We illustrate these in Figure^ The boundary curves of 
each orbit can be worked out analytically. We give the equa- 
tion in Appendix A. For the fast tumbling bar, the orbits 
range from direct to retrograde. As the tumbling frequency 
approaches zero, all orbits become rectangular boxes. By 
properly weighting orbits with different size and axis ratio, 
F66 was able to produce a self-consistent elliptical bar. Its 
distribution function, which is positive definite, has the sim- 
ple form 



/(J) = /o(l-J)-', if 0<J<1, 
— otherwise. 



(7) 
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Figure 3. The different coordinate systems used, x = along 
the heavy dashed Une. 



The analytical expression for J is fairly lengthy, but is es- 
sentially a linear combination of the two energy integrals, 
and is given in equations (43)-(45) of F66. 

The system has a very simple velocity distribution. In 
the frame rotating with the bar, stars stream in the retro- 
grade sense with elliptical streamlines concentric to the bar's 
boundary. The velocity ellipsoid is generally anisotropic with 
the major axis either in the X orY direction. The amplitude 
of the streaming motion is linear with respect to the coordi- 
nates, and the amplitude of the dispersions decreases to zero 
with a radial dependence cr oc (1 — 0)^''^, identical to that of 
the surface density given in eq. (^. A selfconsistent Freeman 
bar model is fully specified by its mass M, the semi-axes a, 
b and the pattern speed SI. One can construct a sequence of 
bar models from needle-shaped to axisymmetric by varying 
the axis ratio b/a. By varying only fl/Qj one can build a 
sequence of models with the same density distribution but 
different velocity anisotropy and streaming motions. 



2.2 Distributions of the transverse velocity and 
the density in the line of sight 

The microlensing rate of the bar is determined by the distri- 
butions of the transverse velocity and the three-dimensional 
volume density in the line of sight path to a source. We 
assume that the observer is sufficiently far away from the 
center of the Freeman bar that all perspective effects can be 
ignored. Even though this makes the analytical studies here 
not rigorously applicable to the Galaxy, where perspective 
effects are a main clue to the size and orientation of the bar, 
it still provides very useful insights. 

The density distribution depends on the thickness of 
the bar, on the orientation of the bar and on the impact 
parameter of the line of sight of the observer. To modify the 
2D Freeman bar to a three-dimensional model, we "puff up" 



the Freeman bar by a small constant thickness A, so that 
the volume density distribution is given by 

E(x,y) 



piX,Y,Z) = 



A < min(a,6). 



(8) 



ZM use a Galacto-centric Cartesian coordinate system 
(a;, y, z) with the x-axis pointing away from the Sun, making 
an angle a with the bar's long (X) axis. Similar but more 
convenient is the following coordinate system (i, y, z) with 



Cy, 



C 



q ^ ) sin a cos a 



q cos^ a + q ^ sin a 



(9) 



so that it is obtained by applying a shear to the original 
{x, y, z) coordinate in the line-of-sight direction, as illus- 
trated in Figure ^ In these coordinates the face-on density 
of the Freeman bar can be rewritten as 



E(i,y) = 



3M 
2Tiah 



-^0 



(10) 



where 2i/o and 2do are the projected width and the central 
depth of the bar. Both quantities are functions of a, b and 
a, given by 



2/0 = (C167) ^ 



with 



7(0,17) 



— 1 2 2 

; q sm a + g cos ol. 



(11) 



(12) 



A line of sight with an impact parameter y intersects with 
the bar's near end at i = —d{y) and far end at £ = +d{y), 
where 



div) = do 1 



y_ 

vl 



(13) 



so that 2d(y) is the depth of the bar at y. For an axisym- 
metric disc, a = b — yo = do — ^ d(y)^ + y^, and 7 = g = 1. 

The transverse velocity distribution depends on the ve- 
locity field of the bar, the orientation and the impact param- 
eter of the line of sight. It is given by a box car distribution: 



/(J)dv 



(14) 



otherwise. 



if \vy ^ v\ < w, 



where the local transverse streaming velocity w is a linear 
function of {x, y), and the half width w is \/3 times the local 
transverse dispersion, which is proportional to (1 — 0)2. v 
can be split into a component v proportional to x and a (less 
useful) component proportional to y. More rigorously 



V 



V ,2; W /—I 



d{y) 



d^y) 



(15) 



t The non-axisymmetric perturbation in the COBE/DIRBE map 
of the Galaxy has a thickness no less than the thin disc (A > 
2 X 200pc) and semi-minor axis no greater than the corotation 
radius (b < a < 4kpc) so the requirement that A/b <JC 1 is barely 
satisfied. 
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where we have rescaled the x coordinate by d{y), the y co- 
ordinate by yo and the velocities by cr(j/), which is defined 



(16) 



a{y) = ao 1 - % 



and (TO — a{y — 0) is the central transverse dispersion. The 
parameter A is defined by 

b r2 1 -iiiiix (17) 



, /b il \ Umax 

\a \lj J (70 



where Umax is the transverse rotation speed at i = ±(io along 
the 1/ = line. So A is a dimensionless global (independent 
of y) indicator of the amount of rotational support in the 
system (similar to the conventional Vrot/(y parameter, e.g., 
Binney & Tremaine 1987). Both oq and A are lengthy func- 
tions of the bar axis ratio and pattern speed, which are given 
in Freeman (1966) in the case that the bar angle a = 0° or 
90°. The dependence on a is then easily worked out with 



(18) 



For our purpose it is only important to know that for a non- 
rotating bar A = and the velocity ellipsoid is aligned with 
the bar; for a Jacobi-type rotator (fi = flj), stars have an 
isotropic dispersion on top of a solid body pattern rotation 
with respect to the rest frame. 

More relevant to the microlensing is the lens-source rel- 
ative transverse speed V and its distribution F(V), which 
is a convolution of the distribution P{vy) of eq. ([14|) for the 
transverse velocity distribution of the lens or the source. 



o-o(a) 
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f cos^ 




1^ sin^ 





FiV) = / P{Vy) [PiVy -V)+ P{Vy + V)] dv y 



(19) 



where we have integrated over the source velocity Vy. The 
integration can be carried out, but the result is somewhat 
lengthy, and is given in eq. (pll). We remark that F{V) 
satisfies the normalization °° F(V)dV = 1. 



2.3 Optical depth 

Following ZM, we assume that the sources are distributed 
with a number density proportional to the mass density ps, 
which is in turn proportional to Es . Then the lensing optical 
depth averaged over all sources along the line of sight with 
impact parameter y is 



\-d(y) 



where 



dxs Es 



+d(y) 
dTs / / dxs Es 



-d(y)<Si<*B 



-d(y) 



dTs 



AtvG 



(xs - xi)pidxi, 



(20) 



(21) 



is the optical depth for a source located at {xs,y,0), and 
Pi is the volume density at the lens position. Here the sub- 
scripts s and ; denote quantities of the source and the lens, 
respectively. For the Freeman bar: 



USGMir_y^^-'^' 
157rc2 A 7 I yi 
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3 EVENT RATE DISTRIBUTION 
3.1 Definitions 

The duration of a microlensing event is related to the lens 
mass m and to the relative distance Xs — xi and velocity V 
between the lens and the source by (cf. eq. 0) 



lAGm 



(23) 



where we have made the approximation D ^ Xs—xi for self- 
lensing of a far-away bar. If one fixes the lens mass m and 
the source position is, then the probability of observing an 
event with any duration is the optical depth dr^ (cf. eq. Gil) . 
But only the fraction F{V)dV of the lenses with relative 
velocity V toV + dV contributes to the event with duration 
t to t -f- dt, where FiV) is the relative velocity distribution 
given in eq. ( |l9| ) . So the contribution dV to the event rate is 



dV= — dTsF{V)dV, 

TVt 



(24) 



where 2/{-Kt) is the average frequency of an event with time 
scale t, and 2/n is the ratio of the Einstein diameter to the 
area in dimensionless units (Paczyhski 1991). It follows that 
the differential duration distribution is given by 

dP 2 In 10 



dlog t 



Tit 



-dTsF{V)V. 



(25) 



Just as for the case of the optical depth (|20[), the observ- 
able rate should be averaged over the source distribution 
along the line of sight. We define the microlensing duration 
distribution profile /(log t) normalized by the optical depth 

as 

(26) 



/(log*) 



rd log t 



so that 



/(logt)dlogt = 



(27) 



where P is the total event rate. Then the source density 
averaged duration profile is given by (cf. eqs |^ and |l|) 

+d(y) X3 

/ dX, / dXi(Xs-5ii)EsEiVF(V) 
/(logt)^^ --^;^^ . (28) 

irt +d{y) Xa 

/ dXa J dxi(Xs— Xl)EsEl 
-d(y) -d(y) 

Note that /(logt) is independent of the bar thickness A. 



3.2 Scaling relations 

The normalized event duration distribution /(logt) should 
have the dimension of a to-be-defined typical frequency u, 
multiplied by a function of the dimensionless time ut. We 
write this concisely as 



/(logt) 
where 



vgx{iyt), 



(To 



'-do)' 



1-y- 

2 

yo 
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Figure 4. Upper panel: the event duration distribution for sta- 
tionary Freeman bars with various axis ratio, angle, and impact 
parameter (fe/a, a, y/yo)', the curves are labeled in the sequence of 
their peak positions (from top left to bottom right). The curves all 
have the same shape, in particular, the same asymptotic power- 
laws for very short /(log t) oc or very long events /(log t) oc t~^. 
Lower panel: a Jacobi-type rotating model (O = S7j) compared 
with a non-rotating one (Q = 0), and the simple double-power- 
law approximation to the distribution (cf. eq. |35| ). The rotating 
models have a narrower distribution and higher rate. 



and 1/:/ is a typical time scale for the events, A is the ratio 
of rotation vs. random motion in the model, and (?a(m) is a 
dimensionless function of fi whose functional form depends 
only on A, which is a dimensionless function of b/a, fi/Slj, 
and Q. The exact derivation of how f(logt) depends on A is 
given in Appendix B. 

Eq. ( [29[ ) can be understood intuitively as follows. Con- 
sider a line of sight with impact parameter y to a station- 
ary Freeman bar model. The lens-source distance {xg — Xi) 
scales with the bar's depth d(y), and the transverse velocity 
V scales with the dispersion cy{y) (cf. eq. |l^). So the Ein- 
stein radius Re oc [md{y)]^^^ and the event rate scale with 

1 _ V ^ a{y) 



, — „ , , ,-r7T OC As a result, the event duration 
distribution depends on the bar density parameters, angle 
and impact parameter only through u. For rotating bars, the 
profiles also depend on the dimensionless quantity A. 

A surprising but direct result of the scaling relation 
eq. (^^ is that all non-rotating Freeman bars have the same 
"shape", as illustrated in Fig. ^. This means that although 



the log t vs. log / diagram of models with the same A have 
different median event duration and total rate, they will 
coincide after shifting the zero points in both axes by a 
constant logi^. In particular, there is a one-to-one corre- 
spondence between a non-rotating Freeman bar with a non- 
rotating Kalnajs (1976) disc, which is an axisymmetric ver- 
sion of the Freeman bar. 

Interestingly the event duration distributions of all 
Freeman bars with the same A also should have the same 
"shape" . The shape is a function of A only. For non-rotating 
models \ = Q. = Q. One can always find a Kalnajs disc with 
the same rotation-to-dispersion ratio (A) and the same shape 
of the microlensing event duration distribution. Just as for 
the dimensionless A (^,«, the width of the event du- 
ration distribution depends on the bar angle and axis ratio 
directly, but depends on the size and mass of the bar only 
through the normalized pattern speed Q/Qj{AI,a,b). 

The shape of the event duration profiles is also invari- 
ant when the same bar model is viewed at different impact 
parameters because the global parameter A is independent 
of the impact parameter y (cf. Fig. 



3.3 Asymptotic behaviour 

Now we consider the asymptotic behavior of the event time 
scale distribution. The asymptotic power-law profiles are ev- 
ident in Figure ^. This is because very short events come 
from a pair of lens and source which are very close in dis- 
tance, and very long events happen if they are very close 
in proper motion. Mao & Paczyriski (1996) show that this 
results in two generic power-law profiles for /(logt) at very 
small or large values of t for three-dimensional models. Sim- 
ilar to these authors, we found the following asymptotic re- 
lations for the two-dimensional Freeman bar, |^ 



/(logt) 



108 In 10 , ^,3 

357r 
108 In 10 



0, 



357r 



oo, 



(31) 
(32) 
(33) 
0.94 



where ^ is a dimensionless function of A with f(0) 
(cf. eq. H). 

The full profile of /(log t) can be constructed, approxi- 
mately, by interpolating between the asymptotic relations 



filogt)^gx{iyt) 



108 In 10 , ,3 



1 



-5p 



(34) 



(35) 



where p is a measure of the width of the distribution. Both 
p and 5 are weakly decreasing functions of A. We find that 

p(A) ^ 0.36 f 1 - ^ V ^(\)^O.m(i + ^) ',(36) 



3 for very long duration, /(logt) generally is proportional to 
(i^t) where dim is the dimension of the system. So for three- 

dimensional bars, ■i/(logi) oc (ut)~'^ when t —> oo. 
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together with eq. (|35[), give a reasonably good approxima- 
tion to /(logi), good within 10% for non-rotating Freeman 
bars = A = 0) (cf. Fig. ^). Using this interpolation, we 
find that logt has a mean at approximately —0.14 — logi^ 
dex with an rms width approximately 0.32 dex; both the 
mean and the width are also weak functions of A. 

The function /(logt) is independent of A or rotation 
for very short events because the lens and source of these 
events are very close in distance, so their relative rotation 
speed is always zero. The short duration events are partic- 
ularly useful for constraining the mass of the lenses in the 
bar because lenses in the foreground disc are generally suf- 
ficiently far away from the sources in the bar that they do 
not contribute significantly to the short events. 

Compared to non-rotating models, the /(logt) profiles 
are slightly narrower for models with increasing rotation. 
This can be understood as a smaller spread in the veloci- 
ties of lens and source in rotating models leads to a smaller 
spread of the event duration. However, the variation of the 
width is small (see Fig ^ , and to good approximation Free- 
man bars all have very similar /(logt) profiles. 



4 THE TOTAL EVENT RATE AS FUNCTION 
OF THE BAR PARAMETERS 

4.1 The size and mass of the bar 

The total event rate per optical depth is given by eqs. (|2 
and (M): 



- = / ugx{i^t)d{logt) = 



(37) 



where is a dimensionless increasing function of A. We find 



Ci(A) ^0.4 1 + 



A^ 



(38) 



Eq. ( p'q) shows that £ is a function of v and A only, and it is 
proportional to v. We now show that several useful relations 
follow from eqs. (^) and (p9|). 

The event rate scales with the lens mass m, the bar mass 
M, size L = y/ab, and thickness A. We define the following 
characteristic lens-source velocity V , Einstein radius R%, 
and frequency v* , 



V = 

Re = 
* 

V = 
Then 

r 

— oc u ■- 

T 



GM 
IT' 

Gm , 



(39) 



R*E 



(40) 



Upon substitution in eq. (^^ , we find the following scaling 
relation for the optical depth and the event rate. 



GM 



— — , FocF* 



G M 





0.0 0.2 0.4 0.6 0.8 1.0 
impoct porometer y/yQ 



Figure 5. The optical depth, event rate, and event rate per opti- 
cal depth (scaled with values at j/ = 0) as functions of the impact 
parameter y. The behavior is the same for all Freeman bars. 



For a compact bar with high M or low ab, the velocity dis- 
persion is high, so the event duration is shortened, and the 
event rate is increased as expected. 



4.2 The impact parameter of the line of sight of 
the observer 

Moving away from the center within a model, the event rate 
drops faster than the optical depth, mostly as a result of a 
lower escape velocity at large radii, so that the mean event 
duration oc u'^ shifts towards larger values (cf. eq. |29| and 
Fig. ^). Specifically: 

3 



p7 



y_ 

yl 



1-^ 
vl 



(42) 



(43) 



(44) 



(41) 



However, the shape of the duration distribution profile is 
always independent of the impact parameter (cf. eq. |2^ and 
Fig I). 



4.3 The bar axis ratio and orientation angle 

For the same bar axis ratio, and the same amount of rota- 
tional support (say A = as for stationary bars), the optical 
depth and the event rate increase if the bar points closer to 
the line of sight, but the events are longer (cf. Fig. ^). This is 
because when the bar points towards us, it is longest in the 
line of sight, and smallest in transverse velocity (motions are 
primarily along the long axis). As a result, both the Einstein 
radius and the event duration are at their maximum. Fig- 
ure ^ also shows that the event rate is significantly higher 
for a rotating model than for a stationary bar with the same 
axis ratio and bar angle. 

At a fixed angle less than 45°, there is an optimal axis 
ratio (near b/a = tana) for the bar to yield the highest rate 
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Figure 6. The optical deptli (upper panel), event rate F (middle 
panel) and event rate per optical depth E (lower panel) as func- 
tions of the axis ratio of the bar at several values of the bar angle 
for a line of sight through the center of the bars. The shaded ar- 
eas show the spread due to different pattern speed. V and r have 
been scaled with F* and r* (cf. eq. 



and optical depth (see also Fig. 2a of ZM, and Fig |6| here). 
At angles larger than 45°, the rate and the optical depth are 
highest for the axisymmetric model. 

Finally, the event rate F largely follows the trend of r 
as expected. This is because 



F (M^ 
— (X u cx 
r 



yl 



(45) 



where 7(0, ^) is given in eq. ([l^), and is generally a "weaker" 
function of the bar parameters (M, a, h) and the line-of-sight 
parameters (a, y) than is the optical depth 



r oc M7 



(46) 



5 TOTAL EVENT RATE FOR MODELS WITH 
THE SAME OPTICAL DEPTH AND 
PROJECTED DENSITY 

Now we return to the question raised in the Introduction: 
how does the event rate vary for models with the same opti- 
cal depth and projected density? This is interesting because 
the projected density from COBE maps and the observed 
optical depth of a bar do not uniquely specify the volume 
density of the bar. Also the velocity structure of the bar 
is not uniquely constrained by self-consistency. If one tries 
to derive the lens mass by fitting the observed microlensing 
event rate with any specific model, the result will be model- 
dependent. Here we try to study this systematic effect. 



5.1 Pattern speed 

First, one can generate a sequence of Freeman bars with the 
same density but different pattern speeds. 

Fixing the bar's volume density and the angle to the 
bar, the event rate is generally (at least for models which 
are bigger in depth than width) larger for rotating models 
than for non-rotating models (cf. Fig. ^ and Fig. g). This 
is because rotation increases the relative lens-source speed. 
The effect can be as strong as changing the event rate by 
25%, as illustrated in Figure ^). 

Since the event distribution is narrower and the rate is 
higher for rotating models (Fig. when fitting these mod- 
els to the observed event distribution, one would derive a 
broader mass distribution with more massive stars as com- 
pared to a slightly narrower spectrum with more low mass 
objects which would be inferred if, instead, one adopts non- 
rotating models (cf. eq. ^ . The change of the mass spectrum 
occurs because the observed spread in log t is partly due to 
the spread in lens mass and partly to the width in /(log t) 
of a single lens mass model. 



5.2 Shear 

One can build a sequence of bar models by shearing along 
the line of sight (Fig. ^). Each model of the sequence has 
the same surface density and optical depth for an observer 
located sufficiently far away from the bar, but all have dif- 
ferent event rates. At least in the limit of a bar sheared to a 
needle-shape, it would have infinitely shallow potential well, 
hence infinitely small lens or source velocity, and would pre- 
dict very few events, but these would last extremely long. 

The sequence of sheared bar models is characterized by 
a ratio between the central depth of the bar and the pro- 
jected size of the bar, do/j/o- The sequence that contains the 



© 1998 RAS, MNRAS 000, 



9 




oxis rotio q = b/o 



0.50 
0.40 

0.30 

0.20 ^ yo/do=0-8>b/a 

V 

0.10 
0.00 



0.2 0.4 0.6 0.8 
oxis rotio q = b/o 



1.0 



0.50 



0.00 




0.2 0.4 0.6 0.8 
oxis rotio q = b/o 



1.0 



Figure 7. The event rate as a function of the axis ratio for a 
sequence of three bars related by a shear along the line of sight 
with a fixed depth to width ratio do/yo = I/7. The shaded areas 
show the spread by bars with different pattern speed. The rates 
are all for a line of sight through the center of the bar. The optical 
depth and the surface density is invariant in the sequence (cf. 
eq. 



Figure 8. A bar model made by shearing an axisymmetric model 
along the line of sight direction. For any line of sight, the shear 
moves the centroid of the mass distribution without changing the 
amount of lens and source, and their relative distance in the line 
of sight, so the projected density and the optical depth do not 
change from model to model. The lower panel draws schematic 
distributions of the transverse speed at the center (the plus sym- 
bol) for a non-rotating bar or disc. The elongated bar has a shal- 
lower potential and lower escape velocity oc y^5bar{0) than the 
disc (oc 2<I>;iisc(0)), where <& is the gravitational potential. 



do /yo > 1 produces more optical depth than an axisymmet- 
ric model, while a sequence with do/yo = I/7 < 1 produces 
less. This is because r/r* oc I/7 = do/yo (cf. eqs |l^and po|) . 
The observed high optical depth towards the bulge suggests 
that do/yo > 1 for the Galactic bar, that is, it is longer in 
depth than across. 

Among bars related by a shear transformation, the rate 
becomes smaller when the bar's axis ratio decreases {b/a 
0), as a result of shallower potential and smaller velocities. 
The rate is largest when the axis ratio is maximum and the 
bar points towards or perpendicular to the observer (Fig 
The velocity in a needle-shaped bar is close to zero, and so is 
the event rate. However, excluding such extremely elongated 
bars which are likely unstable, we find that for bars with 
do/yo = 1/7 > 1 the rate is only a weak function of axis 
ratio; the fractional change is less than 20%. 



axisymmetric model satisfies do/yo = ifj- A sequence with 

^ For these models b/a = tan^ a and the optical depth is the 
same as for the axisymmetric model, while bars with maximum 
optical depth have b/a = tana, and the optical depth is enhanced 
by l/sin2Q (ZM). 



6 THE RATE OF EXTREMELY SHORT 
EVENTS 

As mentioned in §3, the lenses in the bar dominate the events 
at the short duration end. The event distribution is a power- 
law for small t (cf. eq. ^2|). To quantify the rate of these 
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Figure 9. The rate of short events as a function of the axis ratio 
for a sequence of bars related by a shear in the line of sight with 
a fixed depth to width ratio 7=^ = 1; 0-8, 0.5 from the top 
shaded areas to the bottom. The shaded areas show the spread 
by bars with different pattern speed. The rates are all for a line 
of sight through the center of the bar. The optical depth and the 
surface density is invariant in the sequence, (cf. eq. Ed). 



events, we consider the accumulative event rate for events 
shorter than t, 

r{<t) = T I /(logt)d(logt) (47) 



Kt 



t 



0, 



where for the Freeman bar the coeficient K is given by 
36 4 



K 



357r 



(48) 



It follows that r(< t)\t-,o is also a power-law of t with slope 
3 (cf. Mao &Paczyhski 1996). The expression of K shows 
two properties: (a) r(< t)\t^o is independent of the rotation 
parameter A. (b) r(< t)\t^o is very sensitive to u. 

It is straightforward to repeat the analysis of §4 for the 
total rate but now for K, since 



K 4 /r 

Qti V OC I — 

r \T 



1 + 



For example, 



K 4 
V OC 

r 



\ab J 



1 - 2- 

2 

Vo 



(49) 



(50) 



where 7(a, ^) is given in equation (|l^). The general result 
for extremely short events is that both r(< t)|t^o and K are 
very sensitive functions of the bar parameters (M, a, h, Q., a) 
and the impact parameter y. 

The value of K can easily vary by a factor of (1.2 — 
1.3)** ~ 2 — 3 among models with the same projected density 
and optical depth (cf. Fig. This implies that the mass 
moment of the lenses in the bar is poorly determined 

unless kinematic data are used to break the degeneracy of 
models with the same projected density and optical depth. 
To clarify this point, it is helpful to rewrite eq. (kq) as 



/ 1 \ 357r dl r , ,n 



Figure 10. Dependence of lensing properties on the impact pa- 
rameter of the line-of-sight. The horizontal axis is chosen to be 
the run of velocity dispersion from the maximum to zero instead 
of the conventional run of the impact parameter y from the center 
outwards (but cf. eq. and fig. ^) because in this way the curves 
are likely to be valid also for more general bars as well. The event 
time scale (dashed line) increases slowly from the center outward, 
but the rate of very short events (heavy dot dashed line) falls off 
sharply. The optical depth and the total rate are in between these 
two. 



So the extremely sheared and/or rotating models, which 
have a lower ctq, give a higher estimate of (;^). This can 
signficantly influence the determination of the lens mass 
spectrum at the faint end (Mao & Paczyriski 1996). 

r(< t)|t^o follows the general trend of the total rate F 
with the bar mass and size and the impact parameter but 
with a steeper gradient (cf. Fig. ^o|). The dependence on the 
bar angle and the axis ratio is somewhat different from the 
total rate (cf. Fig. ^ and pT[ ). r(< t)\t^o is largest for an 
oblate model (cf. Fig. |ll]) , which does not exactly follow the 
general trend of the optical depth. 

The main difference with the behaviour of F is that 
F(< t)\t-,o decreases with rotation while the total rate F 
increases. This is because as more kinetic energy is put into 
rotation, the dispersion cro is lowered, so if oc oc (cf. 
eq. M) is lowered. 



T CONCLUSIONS AND IMPLICATIONS FOR 
THE GALACTIC BAR 

We have evaluated the event rate, as well as the event du- 
ration distribution and the optical depth, for a family of 
bar models known as the two-dimensional Freeman bars. 
We presented several analytical formulae which show the 
dependence of the optical depth and event rate on the bar 
mass, size, axis ratio, pattern speed, the bar angle and the 
impact parameter with a given line of sight. Models with 
the same optical depth and projected density make slightly 
different predictions on the event rate and the event dura- 
tion distribution for a fixed lens mass. Here we consider the 
implications for the bar in the center of the Galaxy. 



(51) 
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Figure 11. The rate of short events as a function of the axis 
ratio for a sequence of bars with a fixed bar angle. The shaded 
areas show the spread by bars with different pattern speed. The 
Jacobi-type rotators (dashed line) are also isotropic, so K (x 
is the same for the two bar angles. The rates are all for a line of 
sight through the center of the bar. 



7.1 Event rate as functions of the bar parameters 

Can we reliably generalize the results for the event rate of the 
Freeman bars to the Galactic bar? Realistic bars are three 
dimensional and their density structure cannot be modelled 
by the Freeman disc. Furthermore, the orbits in a realistic 
bar potential have more variety (e.g., Pfenniger 1984) than 
in the Freeman bars. For the Galactic bar, it is most inter- 
esting to consider models which can fit both the COBE map 
and the observed optical depth. The Freeman bars certainly 
are too simple to model these. 

Nevertheless the Freeman bars should be able to of- 
fer some insights and serve as a reference to more complex 
systems. The trends found here for the event rate can be 
extrapolated with some modifications. In particular some of 
the results are due purely to geometry, e.g., the event rate 
and duration as functions of the bar angle. Some are generic, 
e.g., the asymptotic power-laws of /(logt); the slope at the 
short duration side is 3 for both two-dimensional and three- 
dimensional models. Others, e.g., the radial increase of the 
event duration, and radial fall-ofi' of the event rate and op- 
tical depth are natural results of a density gradient and the 
radial fall-off of velocity dispersion; the gradient of ^ and 
r will depend sensitively on the density gradient. Because 
the event rate is related to the bar dynamics by several inte- 
grations in the lens-source velocity and volume space, which 
smooth out details of the phase space, only geometrical or 
generic effects will show up in the event rate. So although 
true three-dimensional bars have phase-space distribution 
functions that may be quite different from those of the sim- 
ple Freeman bars, we expect that the latter will match the 
basic microlensing physics in the more realistic models. 



7.2 Uncertainty of predicting the lens mass 

How severe is the non-uniqueness of the lens mass spectrum? 
The common way to derive the mass spectrum in the COBE 
bar using microlensing is to compute first an event distribu- 
tion /(logt) for a lens with mass m with a given dynamical 
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Figure 12. Distribution of the Einstein radius crossing time 
of observed microlensing events towards the bulge. Only data 
from earlier publications of the MACHO (-1- symbols) and OGLE 
(x symbols) teams are shown. Also shown are predictions for a 
line of sight towards the center from some simple models with 
the lens mass being 0.075,0.15 or Q.SMq and the detection ef- 
ficiency crudely taken into account. The three underlying dy- 
namical models are an oblate rotating Kalnajs (1976) disc with 

semi-axes a = b = 1.5kpc ( ), and a simple Freeman (1966) 

bar pointing a = 30° from our line of sight with semi-axes 
(a, b) = (1/ tan a, tan a) X 1.5 = (2.5, 0.8)kpc in the rotating case 

( ) and the non-rotating case (- -). All models have the same 

mass M = 2 X IOIOM0, projected size 10° ~ 1.5kpc and line-of- 
sight depth ~ 1.5kpc, so to yield almost indistinguishable maps 
of the surface density and optical depth. 



model of the bar, and then to convolve it with such a mass 
spectrum that the result fits the observed distribution (e.g., 
Mao & Paczynski 1996). Any non- uniqueness of the dynam- 
ical model propagates to /(logt), which in turn propagates 
to the mass spectrum. 

The current best estimates of the lens masses, which 
use realistic models for the Galactic bulge or bar, indicate 
masses near 0.15Mq, with an uncertain but small < 30% 
fraction below the hydrogen burning limit (Zhao, Spergel 
& Rich 1995; Zhao, Rich & Spergel 1996; Han & Gould 
1996; Han & Lee 1997). If these results were insensitive to 
the adopted density profile and flattening of the bar, then 
a comparable result would be expected from the Freeman 
bars. But if one simply scales up the Freeman bar models 
to the often quoted mass and size of the Galactic bulge/bar, 
and makes a straightforward prediction of the event distri- 
bution (cf. Fig. one finds, quite surprisingly, that the 
typical lens mass is now boosted to around O.8M0, a factor 
of 3 — 5 times larger than the value in more realistic mod- 
els. Fig. ^ also shows that the theoretical predictions can 
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not distinguish very well between a model with O.O75M0 
lenses and a model with the lenses being twice as massive 
when different bulge/bar axis ratio, orientation and rota- 
tion speed are considered. We conclude that the fraction of 
brown dwarfs in the Bulge is also sensitive to uncertainties 
from the dynamical model. 

There are two sources of non-uniqueness of the dynam- 
ical model. The first is the shear transformation. Its efi^ect 
on the distribution f{\ogt) of event durations is generally 
weak, typically at the 10% level in terms of the total rate 
(cf. Fig. 0). Also, as detailed analysis of the COBE map in- 
dicates, the shear transformation generally leaves detectable 
signatures in the left-to-right asymmetry map of the COBE 
map due to perspective efi'ects (Binney et al. 1997, Zhao 
1997 and references therein). The second source of non- 
uniqueness is the variety of velocity structures that are pos- 
sible in the same bar. This non-uniqueness cannot be con- 
strained from the COBE map. As shown also in Fig ^ it can 
create typically about 20% difference in the event rate when 
one compares a stationary bar with a Jacobi-type rotating 
bar. This kind of non-uniqueness is potentially dangerous. 

Most relevant to the lens mass function in the bar is the 
event duration distribution at the short event side. In par- 
ticular, the mass moment (;^) (cf. eq. ^l|) is very sensitive 
to the velocity structure and the shear, and can easily vary 
by a factor of 2 — 3. This translates to an uncertainty of the 
lens mass (particularly at the lower end) by a factor about 
1.4 — 1.7. This implies that any prediction on the faint end 
of the lens mass function in the bar is sensitive to details of 
the adopted dynamical model of the Galactic bar. 

It is possible to reduce the non-uniqueness by fitting the 
stellar radial/proper motion velocity of the bar stars. In fact 
for the Freeman bar, it is possible to lift all the degeneracy 
if one can measure the line of sight mean streaming motion 
and velocity dispersion. For more realistic three-dimensional 
models, the non-uniqueness of the phase space structure is 
more complicated, and must be constrained by detailed nu- 
merical modeling. Generally one would need to fit radial ve- 
locities and proper motions at a number of positions of the 
bar. Combined with the COBE map and the microlensing 
optical depth, these data will strongly if not uniquely limit 
the axis ratio, the angle of the bar and the pattern speed. 
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records the hospitality of the Max-Planck Institute in Garch- 
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APPENDIX A: ORBIT BOUNDARIES 

In the frame corotating with the bar the curves bounding 
the area filled by an orbit of the Freeman bar can be written 
in parametrized form using eqs (31-33) of F66. We find 



X(<^) 



y(0) 



1 ± 



1 ± 





kiv^Y 










+ V 


1 - 

2 









(Al) 



(A2) 



where the ± sign indicates the outer and inner boundaries, 

fj. — Ac cos (j), V = kaAaSilKj), (A3) 

and the parameter is an angle ranging between and 27r. 
All quantities have the same meaning as in F66. 



APPENDIX B: DERIVATION OF THE EVENT 
DURATION DISTRIBUTION 

The distribution of the relative lens-source speed is inte- 
grated to be (cf. eq. jl9{) 



F{V) = max[0, QiV)] + max[0, Qi-V)], 
where 



(Bl) 



'lWsWiQ{V)= min [ws,wi + vi^ - V] (B2) 
+ mm[ws,wi — vts + V], 

and Vis = vi'-Vs- Apply the following change of variables to 
eq. M) 



(B3) 
(B4) 



Xs = d{y) smiles, xi — d{y) svaipi, 
define A and v as in eq. (^ij), and write 
5 = sin-!/);, — sin?/);, n = ut, 
then 

Ws,i = VSo-{y) COS ip 3, 1, Vs,i = Xa{y)sm^s,i, (B5) 
and 

7o' 



2^j.-^Vd. (B6) 



It follows that (cf. eq. ||) 



where 
7i(A,m) 

and 



(B7) 



/2 pi>„ 

dtps / dtpi 5 cos^ ipa cos^ '(/'I (B8) 

r/2 J-tt/2 

X (max[0, Fo(A, fi, tps,il'i)] 

-|-max[0,Fo(A, -^i,ips,ipi)]), 



r''^ f^" 2 2 32 

I2 = dipB / drpi S cos Tpe cos V"! = 77 ■ (B9) 

Here we have defined 



Fo = {ut)VQ{V) 



(BIO) 



6 cos ips cos ipi 

X (min[cos tpa, cos ipi — _ffo] 

+ min[cos7/)s,cosi/;; -t-i^o]), 



where 



1 / , <- 2v(5, 

Ho= — A5+ 

V3 M 



(BU) 



Clearly 1\ is generally a dimensionless function of A and vt 
only, so that 
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-f{\ogt)=g,{vt) = ^^h{\,vt). (B12) 

The asymptotic relations (cf. eqs ^ and ^) can be derived 
by letting 5 — > and ^ for t — *■ 0, and /i — > oo for 
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